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DISCUSSION OF PLATES WITH BOUNDARY CONDITIONS 
OF ELASTIC SUPPORT 
PROCEEDINGS-SEPARATE NO. 199 


M. Z. v. Krzywosiocxr.—An interesting approach has been presented by 
18 Prof. of Gas Dynamics and Theoretical Aerodynamics, Univ. of Illinois, Urbana, II. 

the author to the problem encountered in the investigation of elastically sup- 

ported plates. It appears that Mr. Fuchs has underestimated the method 

offered by Mr. Ritz. Although this method does not usually bear Mr. Ritz’s 


‘Uber eine neue Methode zur Lésung gewisser Variationsprobleme der mathematischen Physik,” by 
W. Ritz, Journal fiir reine und angewandte Mathematik, Vol. 135, 1909, pp. 1-61. 


name, the method has furnished many successful solutions in the theory of 
plates. 

In 1908 Mr. Ritz proposed an approximation method for solving the 
problems of a variational character.“ He proved that under some conditions 


4 “Variational Methods for the Solution of Problems of Equilibrium and Vibrations,” by R. Courant, 
Bulletin, Am. Mathematical Soc., Vol. 49, 1943, pp. 1-23. 


the minjmizing sequence converges (but stated that the limit of convergence 
can be quite different from the sought for solution), and he applied his method 
to plate problems, vibrating-string problems, and the Dirichlet problem. 
Subsequent investigators have made use of Mr. Ritz’s work but have failed to 
make the proper references to the initial investigation.’*!*'!7 A comparison 


146“*The Approximate Solution of Two-Dimensional Problems in Elasticity," by 8. Timoshenko, 
Philosophical Magazine, Vol. 47, Series 6, 1924, pp. 1095-1104. 

16 “Strength of Materials,’’ by 8. Timoshenko, D. Van Nostrand Co., Inc., New York, N. Y., 1930, 
Vol. 1, p. 326, Paragraph 7.1. 

17 “Theory of Elasticity,’’ by 8. Timoshenko, McGraw-Hill Book Co., Inc., New York, N. Y., 1934, p. 
150, Paragraph 43. 


of Mr. Ritz’s work with the subsequent investigations has been made by the 
writer.'* It has been shown that the principle of least work (based on the 


18 “On the So-called Principle of Least Method,” by M. Z. v. Krzywoblocki, Proceedings, First Mid- 
western Conference on Solid Mechanics, Univ. of Illinois, Urbana, Ill., April, 1953. 


method proposed by Mr. Ritz), depending on the choice of the minimizing 
sequence, may yield a rapid convergence (as shown in some of the examples of 
S. Timoshenko) or nonsensical results (as shown by the writer). It can be 
assumed that each particular type of boundary-value problem is character- 
ized by its own most proper minimizing sequence. 

Other approaches (developed by 8S. Bergman) to the problem under in- 
vestigation are the integral-operator method"® and the method of solving the 


1 “Uber Schubknickung von isotropen und anisotropen Platten,” by 8. Bergman, Proceedings, Third 
International Cong. of Applied Mechanics, Stockholm, Sweden, 1930, pp. 82-87. 


boundary-value problem.” Mr. Fuchs did not mention this interesting 


2% ‘Neuere Probleme aus der Flugzeugstatik. Uber die Knickung von rechteckigen Platten bei Schub- 
beanspuchung,” by 8. Bergman, Prikladnaja Matematika i Mekhanika, Vol. 2, 1935, pp. 207-224. 


method, proposed by Mr. Bergman in 1930 and appiied by the writer in 1948.7! 


21 ‘A General Approximation Method in the Theory of Plates of Small Deflection,” by M. Z. v. Krzy- 
woblocki, Quarterly of Applied Mathematics, April, 1948, pp. 31-52. 


Other methods which can be successfully applied are those that were developed 
by A. and R. V. Southwell.**?7 


22 ‘Etude des spectres des équations aux dérivées partielles de la théorie des plaques élastiques,”” by A. 
Weinstein, Mémorial des sciences mathématiques, No. 88, 1937. 
23 ‘On the Unified Theory of Eingenvalues of Plates and Membranes,” by N. Arojszajn and A. Wein- 
stein, American Journal of Mathematics, Vol. 64, 1942, RP. 623-643. 
2% “On the Bending of a Clamped Plate,’’ by A. Weinstein and D. H. Rock, Quarterly of Applied 
Mathematics, Vol. 2, 1944, pp. 262-266. 


25‘Qn a Boundary Value Problem for a Clamped Plate,” by A. Weinstein and J. A. Jenkins, Trans- 
actions, Royal Soc. of Canada, Section 3, 1946, pp. 59-67. 

2% ‘Relaxation Methods Applied to Engineering Problems,” by L. Fox and R. V. Southwell, Trane- 
actions, Royal Soc. of London, Series A, Vol. 239, 1945, pp. 419-460. , ; 

27 “Relaxation Methods Applied to Engineering Problems,” by D. G. Christopherson, L. Fox, F. 8. 


Shaw, and R. V. Southwell, ibid., pp. 461-487. 
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R. H. Wood.?*—An analytical solution which rarely has been attempted 
28 Prin. Scientific Officer, Building Research Station, Garston, Watford, Hertfordshire, England. 


because of its complexity has been obtained by the author. The problem 
involved is the determination of the behavior of rectangular and square plates 
supported on elastic beams. It is not generally understood, however, what 
significance the solution has in relation to structural-engineering research. 

A great deal of attention has been given to the design of multi-story building 


frames both in the elastic field® and by “plastic” methods to determine the 


2 ‘An Economical Design of Rigid Steel Frames for Multi-Storey Buildings,’’ by R. H. Wood, Research 
Paper No. 10, National Building Studies, Her Majesty's Stationery Office, London, Pngland, 1951. 


point of collapse. The next logical step would be to determine what addi- 


%*“\ Review of Recent Investigations into the Behavior of Steel Frames in the Plastic Range,”’ by 
J. F. Baker, Journal, Inst. of C.E., January, 1949. 


tional effeets—in terms of both the behavior at working conditions and the 
enhancement of the collapse load—can be produced by the presence of cladding. 
When there are reinforced concrete floors in a building, it has been found that 
the reduction of stress in the supporting beams can be so large that the resulting 
frame action is considerably modified. In Fig. 10 there is shown a floor being 
subjected to a distributed load of 12 tons. The maximum recorded stress was 
approximately } ton per sq in. The corresponding stress computed for bare 
frame action would be many times this amount. Consequently, the Building 
Research Station at Garston, Watford, Hertfordshire, England, has commenced 
a systematic study of “composite action” by analysis, controlled tests in the 
laboratory, and field tests. It will be shown that the problem which Mr. Fuchs 
has investigated is one of considerable practical importance. 

This particular floor-beam problem has been studied by the writer using the 
method of finite differences (for elastic behavior) and the “limit-analysis”’ 


method (for ultimate collapse), and a general report on the relationship of the 
problem to the study of composite action in multi-story buildings has been 


prepared.** In Fig. 11 there is shown the distribution of vertical reaction 


“Studies in Composite Construction, Part II,"" by R. H. Wood, Research Paper No. 24. National 
Building Studies, Her Majesty's Stationery Office, London, England (publication pending). 


between floor and beam for a square slab with beams of equal stiffness. Fig. 11 


Fic. 10.—Fioor Unper Loap 
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Fic. 11.—DistrisvtTion oF Loap INTENSITY ON THE SupPorTING Beams (yu = 0) 
For Unrrorm Loapine on Square SLAB 


Fig. 12.—Distrisvtion or Moment at Center Section or Stas (u = 0) 
ror Unirorm Loapine on Sqvare SLAB 


has been drawn for the case in which Poisson’s ratio is equal to zero. The 


British Code of Practice® allows the computations for elastic design to proceed 
32 ‘‘The Structural Use of Normal Reinforced Concrete in Buildings,’’ British Standard Code of Practice 
C.P. 114, British Standards Inst., London, England, 1948. 


on a basis of Poisson’s ratio equal to zero, but as an alternative the beam loads, 
as shown by the dashed lines, can be used. It can be seen that Fig. 11 is of the 
same form as Fig. 7 after allowance has been made for the effect of Poisson’s 
ratio. Similar agreement is found in the case of the bending moments across 
the center section (Fig. 12). The parameter relating beam stiffness (£ J) 


to the flexural rigidity of the half slab (v4) is 


It is apparent that, when vy is equal to unity, a ‘“‘twistless’’ case ensues 
(when y is equal to zero)—each strip acting independently and resulting in a 


uniform reaction on each beam of BL per unit length. It is interesting to note 


that the author has found (Fig. 7) that a uniform load on the beams can be 
produced in certain cases even when Poisson’s ratio is not zero—that is, when 
the condition is not that of zero twist. Of even greater interest is the fact that 
the writer has discovered that a uniform beam reaction can ensue at the point 
of incipient plastic collapse when there is a considerable amount of twist 
present. This phenomenon might be termed “reaction-redistribution” and is 
of considerable importance in determining the behavior of composite systems 
near collapse. 

There are some interesting features in the polynomial expression for the 
deflection of a plate. The writer has investigated the deflection of a composite 
slab-beam system and has found that the deflection can be expressed as 


te 
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which is a twistless solution (u = 0), whereas Mr. Fuchs has found, for the 
more general case, 


In order that the support moments shall be zero, Eq. 38 leads to the re- 
striction, 
f 6 yu e= 


and entails the necessary symmetry, 


which leads to the singular solution, 


NZ 


for a rectangle of sides L, and Ly. 

It is interesting to note, however, that, if u is small enough to be neglected, 
there results a greater generality in the possible variations in beam stiffness. 
Commencing with Eq. 37, the bi-harmonic equation for the plate results in 


in which ¢c is not equal to e. Also, since y is equal to zero, 


b 
and by symmetry 
d 


In the twistless case the beam-to-slab reaction is 


dy 


Along x = L,/2, Eq. 43a can be written as 
12 Le c N 


| 
: 
2 
If M, is equal to zero along the side x = + L,/2, 
3 i 
530-4 


so that 


which leads to 


For Eqs. 45 to be applicable simultaneously, it is only necessary to have 
the less restrictive condition, 


It will be noticed from Eqs. 46 that there is now no restriction on the 
magnitudes of the beam stiffnesses. However, in this twistless case if one 
set of beams is strengthened the other set must be weakened. The intensity 
of load on the beams is 


Along z = V. 


Vy 


Along y = 


Consequently, if Poisson’s ratio is small, the reactions will be given ap- 
proximately by Eqs. 47. However, because there is no restriction to the 
singular solution (Eq. 39c) it is easier to interpret the behavior of floor-beam 
systems. 

This conception is enhanced when one attempts to determine if a twistless 
case can exist in continuous slab-beam systems, restricting the investigation 
to continuous rectangular panels supporting continuous uniform loading. 

Assuming 

N (1 — 


at all points, then dw/dy must be a function of y only. It follows that 


The load on the beam, however, is 
EI,24e (44a) 
2%eHI,=—V,= +12L,cN............... (44) 
wet .. (45a) 
Fl. "yc 
and 
(456) 
(46a) 
or 
i+— 
Yz 
1+ — 

Yu 

ver 
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Because this attempted solution represents, in effect, the behavior of 
independent strips, the load carried in the z-direction is not necessarily the 
same ; therefore, 


from which 


= oy! (x) +O0+f' (y) = 


If the load is shared in the proportion, 


= 


Ow _ 


retaining only odd powers of x because of symmetry. 
Since w is equal to zero at the corners, 


48 384 N 48 384 


Along the edge, y = L,/2, the combined reaction per unit length from the 
slabs on both sides is found to be (1 — a) p Ls, and along the edge, x = L,/2, 
the reaction is equal to a p L, both of which are constant. 

Because the beam is virtually built-in (because of continuity), the deflection 
of the long beams carrying this uniform load would be 


L’, LA, 
—a)ply (z- 48 384 


but the slab deflection along y = L,/2 is given by 


(# as + \ 353 Too + 384 


zero 
Eqs. 53 and 54 are identical if 


in which the subscript ¢ denotes continuity, and each beam is now associated 
with a total width of slab of Ly». 
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Similarly, 


from which 


which is the condition relating slab and beam stiffnesses if the twistless case 
is to occur. This condition is valid for all values of Poisson’s ratio, whereas 
the condition was only valid for u equal to zero in the case of isolated panels. 
Of greater interest is the fact that Eq. 56) is exactly the same as the “‘less 
restricted” condition given by Eq. 466. These twistless cases therefore serve 
as a landmark in a general study of the subject. 

Fig. 13 shows the results obtained for the beam loads which appertain to 
continuous, square panels. Fig. 13 emphasizes the importance of studying the 
load distribution on the beams. The finite-difference results contain small 
errors which could be improved by using a finer mesh in the computations. 
The twistless case shows up prominently. 

A brief introduction to the collapse of such slab-beam systems is of interest. 
Fig. 14 shows the collapse of a symmetrical square panel. The collapse mode 
is identified by the fracture lines (lines of full plastic moment) in the slab 
approximately along each center line and by the plastic hinges in the beams 
to form a mechanism. Above a certain critical ratio of beam-to-slab carrying 
capacity, theory indicates that the rectangular mode of collapse is abandoned 
in favor of a diagonal collapse of the slab. The supporting frame does not 
take part in this collapse. There is clearly a limiting load distribution col- 
lected by the beams, which is eventually insufficient to cause the beams to 
collapse. That is, this limiting load distribution is decided by the slab and 
beams working in unison. In general, it cannot be decided intuitively in 
advance and then coupled with a load factor. 

Design on such a basis may, however, be intuitively correct in special 
cases." Consequently, Mr. Fuchs’ analytical solution is a weleome and valu- 
able addition to the information which is required in the study of composite 
structures in general. 
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8. J. Fucus,» M. ASCE.—Interesting comments concerning developments 
% Chf. Stress Analyst, Loewy Construction Co., New York, N. Y. 

in various aspects of the investigation of elastically supported plates have been 
offered by Messrs. Krzywoblocki and Wood. 

Mr. Krzywoblocki’s justifiable criticism is that the writer’s brief comment 
relative to the method offered by Mr. Ritz appears to underestimate the 
method. Under the heading, ‘Introduction: Historical Note,” it was stated 
that the usual application of the Rayleigh-Ritz method involves a minimizing 
sequence 


Wa = CWi + Coe + 


in which each coordinate function y; satisfies the boundary conditions. How- 
ever, functions are not readily available which satisfy the boundary conditions 
of the typical plate, simply and flexurally supported, which has been analyzed 
in detail. This problem can be investigated by use of the strain-energy method 
with dependence on natural boundary conditions; that is, the integral ex- 
pression for strain energy D, (w,, A, x) —with dependence on the parameters 
\ and x—can be minimized with respect to the coefficients c;. This minimum 
evaluation of D, may not be affected intimately by the choice of coordinate 
functions y;, and similarly the evaluation of a scalar such as a characteristic 
value (also derived by integration) may likewise be unaffected directly.- The 
evaluation of other properties, however, such as shear (which depends on third- 
order derivatives) is comparatively directly affected by the choice of functions 
¥;. By the formulation of the problem without the aid of functions which 
satisfy the conditions as imposed, the error involved in the evaluation of a 
property is probably increased. Furthermore, the strain-energy method as 
applied to the case of the typical, flexurally supported plate requires the solu- 
tion of simultaneous linear equations—as is also necessary in the method 
offered by the writer. 

In 1953 a strain-energy analysis of the clamped, flexurally supported plate 
was reported™ by J. G. Sutherland, J. M. ASCE, L. E. Goodman, and N. M. 


“Analysis of Plates Continuous over Flexible Beams,”’ by J. G. Sutherland, L. E. Goodman, and 
N. M. Newmark, Technical Report, Project NR-064-183, Office of Naval Research, Washington, D. C., 


January, 1953. 
Newmark, M. ASCE. Their problem admitted the use of coordinate functions, 
referred to previously, which facilitated the solution by satisfying the imposed 
boundary conditions. Certain functions, formed of combinations of Legendre 
polynomials, were utilized so as to satisfy explicitly the condition of zero normal 
slope around the boundary. For the sake of a comparison of results obtained 
by the alternative methods, some typical maximum moment coefficients derived 
by the strain-energy method* are listed in Table 3, together with corresponding 
coefficients obtained by the use of the writer’s method. No shear coefficients 
were listed by Messrs. Sutherland, Goodman, and Newmark.* The complete 
analysis and numerical evaluation (by the writer’s method) of the clamped, 
flexurally supported plate may be published subsequently. 

The historical note was limited to a few standard solutions that have been 
applied to plate problems. A more extended survey would have required 
mention of the important methods developed by Mr. Bergman. The latter’s 
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(57) 


integral-operator method, alluded to by Mr. Krzywoblocki, provides a solu- 
tion for an assigned, linear, partial, differential equation as a transform of 
arbitrary functions. The method, however, appears to offer no particular 
facility in dealing with the boundary conditions typical of this subject. 
Futhermore, the method requires involved computations.’ The writer’s 
method is considered to be sufficiently general for many applications in engi- 
neering and was shown to reduce to standard solutions for particular cases. 
Mr. Wood’s informative account of research on composite panel action 
conducted at the Building Research Station at Garston (England) indicates 
the alertness of investigators to the implications, in practice, of developments 
in the theory of transversely loaded panels. The coefficients for moment 


TABLE 3.—MaxiIMUM—MOoMENT COEFFICIENTS FOR CLAMPED PLATES? 


VALUE? OF AND A/D 


Coefficient 
0.5 | 0.25 | 0.5 


0.25 | 0.1 


VALUE oF b/a 
By S. J. Fuchs’ method 


1.0 0.75 


—Cmz| 2-0/2, =b/2 1/24 0.069043 0.035715 0.057866 0.097062 


—Cmy| r—a/2, y=b/2 1/24 0.069943 0.026785 0.046045 0.082592 


VALUE oF b/a 


By strain-energy method¢ 


1.0 0.80 


2-0/2, y=b/2 1/24 0.06904 0.03704 0.06010 | 0.09951 


—Ciny| 2-0/2, 1/24 0.06904 0.02963 0.05037 | 0.08775 


¢ At the corners of flexurally supported clamped rectangular plates (u =0) under umform load, p. 
Thus, Mz = Cnz pat, and My = Cm,pa*. Values of parameter refer to single panels, or one half of 
‘between adjacent panels. ¢ As reported by Messrs. Sutherland, Goodman, 


the total ridigity of the beam 
and Newmark. 


and shear which Mr. Wood derived by another alternative method—that of 
difference equations—are of great value in providing the mutual corroboration 
required by the engineer. 

The exact solution for the clamped plate classification can be derived con- 
veniently by recognizing that the normal slope, hence the mixed derivatives 
Wry, Vanish at each point on the boundary. Thus, it is sufficient to set f equal 
to zero in Eq. 38, and Eqs. 43 to 46 follow as shown. 

It is agreed that elastic stress, computed (for example) by the coefficients 
from Table 3, are also of value in consideration of plastic stress. Thus, the 
residual stress caused by plastic action can be computed as the difference be- 
tween plastic and elastic stress at a generic point provided that the residual 
effective shear does not exceed the yield strength in unloading. Thereafter, 
cycles of loading and unloading may be all-elastic; in these cases, elastic 
coefficients are necessary. 


530-9 


| 
—_ 


Flexurally supported panels are important in structural design. It has 
been usual to consider panels to be supported on beams of implied complete 
rigidity (A— ©). The analyses and numerical data which have been made 
available for flexurally supported plates give the designer greater latitude 
and enable him to use supporting beams of arbitrary rigidity. 
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DISCUSSION OF APPLICATION OF THE RELAXATION 
TECHNIQUE IN FLUID MECHANICS 
PROCEEDINGS—SE PARATE 223 


JOHN S. McNOWN', M. ASCE, EN-YUN HSU”, A.M. ASCE, and CHIA- 
SHUN YIH*.—In the several discussions submitted, contributions have 
been made which range from considerations of the philosophy underlying 
numerical methods through matters of technique to comparisons of nu- 
merical and observed results. The fact that only two discussers de- 
scribed the conduct of additional calculations is probably an indication 
that the method has not yet found general use. If the paper and the dis- 
cussions serve to stimulate such use, their primary objective will have 
been achieved. 

Mr. Blaisdell considered primarily the accuracy, or rather the con- 
sistency of the numerical results for flow over a weir. The obvious 
lack of system at one point in his table led to the finding of a typograph- 
ical error in the published results. For the upper nappe, the value of 
y/H for x/H = } should have been 0.442 instead of 0.422. Aside from 
this erroneous value the results of his equation and of the relaxation are 
remarkably close together. 

Althougt Mr. Blaisdell characterizeii a part of the correspondence 
as erratic, Messrs. Paynter and Scott apparently found the general 
agreement between the results observed at the Bureau of Reclamation 
and those of the relaxation almost too good to believe. Actually, these 
data came to the writers’ attention well after both the calculation and 
the manuscript had been completed. Had they been used in the calcula- 
tion, reference would, of course, have been made to them. Mr. Frank 
B. Campbell of the Waterways Experiment Station informed the writers 
of their existence in a letter, and Mr. J. N. Bradley of the Bureau of 
Reclamation submitted a comparison like that of Messrs. Paynter and 
Scott. Unfortunately, his comments were received too late to be pub- 
lished as a discussion. 

A comparison between second differences can be misleading, par- 
ticularly if one set is calculated from a quadratic equation for which the 
second differences are inherently constant. As their order of magni- 
tude in this calculation is approximately H/100, presentation of the sec- 
ond differences as 3- and 4-digit numbers by Blaisdell is misleading. 
Even for his empirical equation, the rounding-off process results ina 
sizable change. 

T. Prof., Univ. of Michigan, Dept. of Eng. Mechanics, Ann Arbor, Mich.; 
Formerly Associate Director, Iowa Inst. of Hydraulic Research, State 
Univ. of Iowa, Iowa City, Iowa. 

2. Physicist, Hydromechanics Lab., David Taylor Model Basin, Navy 
Dept., Washington, D.C. 

3. Research Engr., Iowa Inst. of Hydraulic Research, State Univ. of 
Iowa, Iowa City, Iowa. 530-11 
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The trend of the values for the thickness of the jet is more signifi- ' 

cant. However, the value of 0.5 for x/ H is not large enough for the in- ) 

ternal pressure to be ignored, as is clearly indicated in Fig. 13 of the 

original paper. For x/H~1, the internal pressure head is still more 

than H/10. Thus the conditions of free fall assumed by Blaisdell are 

only approximately valid. The fact that the observed thickness differs 

from the mean of the calculated thicknesses by less than one percent 

seems quite remarkable. 
Messrs. Paynter and Scott objected to the absences of both a broader 

statement of the underlying philosophy of numerical methods and fuller 

numerical detail. In their comment, “Rather than merely being tech- 

niques of numerical integration, ...” the word “merely” is the dis- 

cussers’ and not the writers’. Their more general interpretation of the 

method is in accord with the writers’ views. It had already been amply 

stated in such references as those supplied by the discussers. Omission 

of the relative width (b = 2h/3) of the permeable wall in Fig. 17 was an 

oversight. 

Numerical details were included here and there in the paper more as 


illustrations of what can be done than as exhaustive solutions to practi- 
cal problems. In this connection the computations for weir flow (Fig. 13) 
and for seepage (Fig. 17) are hardly comparable in significance. The 
former is a classic problem of considerable practical importance; the 
latter is an isolated example which is only illustrative of an infinite 
variety of boundary configurations. Furthermore, difficulties inherent 
in the presentation of the numerical results will be obvious to the read- 
er. In the cavity flow (Fig. 14), for example, over 1300 values of the 
stream function were determined. To have gone very far in either of the 
somewhat opposing directions proposed by Paynter and Scott would have : 
unduly extended an already long paper. } 
Boreli and Paynter and Scott described methods and conditions for 
which a mesh other than square is used. Some of these were also dis- if 
cussed in Southwell’s original work on the subject. Aside from the ! 


square mesh, however, only the rectangular which Boreli used to advan- 
tage in calculating seepage toward wells is likely to be applicable to 
problems of fluid flow. The rather distantly related Monte Carlo meth- 
od is interesting and perhaps useful for isolated applications. 

Mr. Boreli’s contention that the formula obtained by Taylor’s expan- 
sion for the residue is often preferable to the one presented in the paper 
is readily conceded. In fact, in using the relaxation technique to deter- | 
mine the cavity shape (Fig. 14), one of the writers used the formula 
from Taylor’s expansion. It must be noted, however, that if the second 
and higher derivatives of the unknown function are large, and the net is 
comparatively coarse, the superiority of Eq. (2) to Eq. (1) in Mr. 

Boreli’s discussion cannot be established, since Eq. (2) is then only a 
rough approximation. All one can say is that as the net becomes finer, 
the accuracy of Eq. (2) is increased more rapidly than that of Eq. (1). 

Mr. Boreli proposed a treatment of singular points by properly 
weighting the near-by legs. Unfortunately, he did not mention how the 
weighting is to be determined. The treatment presented by Mr. Sarpkaya ’ 
for the same purpose is simple and systematic and is a welcome contri- i} 
bution to this troublesome aspect of the relaxation technique. 
530-12 


Although the writers stated that application of the relaxation method 
is feasible only for two-dimensional or axisymmetric flows, recent 
studies*»® have included analyses of problems for which variations in all 
three dimensions are significant. Mr. Roy obtained a rough approxima- 
tion for flow through a transition from a reservoir to a rectangular con- 
duit, and Messrs. Allen and Dennis demonstrated a useful technique 
which is particularly well suited to problems for which the third dimen- 
sion is small compared to the others. A solution of acceptable accura- 
cy for a three-dimensional flow with a free surface, as an extreme ex- 
ample, remains an exorbitant undertaking. 

The discussers indicated that they share the writers’ convictions as 
to the applicability and versatility of the relaxation technique. Mr. 
Sarpkaya’s application of the method to problems of wave motion is ad- 
ditional evidence of its adaptability to a wide variety of potential flows. 
By their contributions, the discussers have expanded and reinforced 
various phases of the study. 


4. Allen, D. N. de G. and Dennis S. C. R., “The application of relaxation 
methods to the solution of differential equations in three dimensions” . 
I. Boundary value problems. Quart. Sr. of Mech. and Appl. Math., Vol. 
IV Part 2, June 1951, pp 199-208. 


5. Roy, S. K., “Development of three dimensional relaxation processes 
for a systematic approach to flow analysis in hydraulic structures” (In 
three parts), Irrigation and Power (Indian Jr.) Vol. 9, No. 1, Jan. 1952. 
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DISCUSSION OF NUMERICAL ANALYSIS 
OF FRAMES WITH CURVED GIRDERS 
PROCEEDINGS—SEPARATE 250 


JAMES MICHALOS', M. ASCE.—Mr. Conwell’s discussion has added 
to the value of the paper, and his very kind comments are greatly ap- 
preciated. 

An interchange of subscripts within the parentheses in the notation 
for carry-over may be desirable. The arrangement chosen seemed con- 
venient, but this is largely a matter of preference. However, because 
this paper is related to another”, and on occasion use of both may be 
necessary for an analysis, some additional remarks regarding notation 
are pertinent. In the latter paper, written and presented previous to the 
paper now under discussion, a prime mark is used to indicate carry- 
over to the opposite end of a member. In the paper now under discussion 
it was necessary to use a more specific designation, such as A-B and 
B-A, for carry-over factors. 

Tables or curves of properties are not at present available for other 
than members of circular-arc plan and constant cross section but, as 
stated in the paper, values of stiffnesses, carry-over facto:’s and fixed- 
end moments can be determined for any shape and variation in cross 
section by application of the shear and torsion analogy presented else- 
where*. In that paper it is shown that moments resulting from loac ap- 
plied normal to the plane of a structure fixed at its ends can be com- 
puted directly by a numerical procedure in a manner analogous to the 
computation of stresses due to direct shear and torsion in a strut. The 
procedure is illustrated in Figs. 6 and 7 of that paper. The reader 
should observe the difference in sign convention used in that paper and 
this one, and should also observe that the y-axis of that paper is the 
z-axis of this one. With these differences accounted for, the fixed-end 
moments in Fig. 7 of the previously published paper become m, = - 8.79 
and mz = -15.01 at A, and mg = - 2.46 and mz = +4.81 at B. 

When account is taken of the differences in axes and signs, the stiff- 
ness at an end A of a member AB, where end A is the right-hand end, 
and the corresponding carry-over factors (in terms of the stiffness) can 
be computed as follows: 

1. Prof. of Structural Eng. and Chairman, Dept. of Civ. Eng., New York 
Univ., New York, N. Y.; formerly Prof. of Civ. Eng., lowa State College, 
Ames, Iowa. 

2. “Numerical Analysis of Continuous Frames in Space”, by James 
Michalos, Proceedings—Separate No. 261, ASCE, September, 1953. 

3. “Laterally Loaded Plane Structures and Structures Curved in Space”, 
by Frank Baron and James Michalos, Transactions, ASCE, Vol. 117, 
1952, p. 279. 
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(Xe), Zc) 

Kra(Cxz)a-a = - + 

(2) 

(xq) 

Kza(Cxz)a-p = x 

(x,), (x,) 

1 c c 


A, yy 


To obtain Kzp and the corresponding carry-over factors, interchange 
all A and B subscripts in the preceding expressions. To obtain and 
the corresponding carry-over factors substitute x for z and z for x 
the expressions. To obtain K,p interchange the A and B subscripts as 
well as z and x. 

If symmetry exists, as in Fig. 7 of the paper* to which reference was 
made previously, or if all elements of the member are parallel to one or 
the other of the reference axes, as in Fig. 6 of the same paper, then both 
P and P; are equal to zero, and A, and A, are equal, respectively, to 
Ax and A,. As an example of the use of Eqs. 1, a stiffness and the cor- 
responding carry-over factors (in terms of stiffness) will now be ob- 
tained for the parabolic arch of variable section shown in Fig, 7. 

First, however, it is necessary that ends A and B of Fig. 7 be inter- 
changed to conform to the nomenclature used in the paper now under 
discussion. Then, for end B (end A in Fig. 7), 


o 2 

= + = 0.0521 EI 
0 (-50)(- 9.42) 

Kp (Cxz)p-_B = 70.04 * $3,086 = 0.0076 EI 
(- 50)(- 9.42 

Krp(Cxz)p_a = -0.0076 EI 
(-50)(50)_ 0285 EI 


where I is the moment of inertia about the z-axis at the crown. 

Tables of fixed-end moments, stiffnesses and carry-over factors for 
any shape can be prepared by writing general algebraic expressions for 
the various terms of the shear and torsion analogy. The tables for mem- 
bers of circular arc plan were prepared in that manner. 

Mr. Conwell outlines a straight-forward procedure for obtaining 
fixed-end moments due to a uniform load applied over a portion of a 
circular girder. Such values can also be obtained directly, as implied 
by Mr. Conwell, from influence lines. Influence lines for fixed-end 
bending moments are shown in Fig. 5. These were drawn by plotting 
values from Table 3. Influence lines for fixed-end twisting moments 
can be drawn in the same manner. To illustrate the method of obtaining 
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fixed-end moments due to partial loading, assume a circular girder AB 
whose central angle, 8, is equal to 90 degrees. Call its arc length L. 
Furthermore, consider this girder subjected to a uniformly distributed 
load, w per unit length of arc, over two-thirds of its length starting at 
end A. The fixed-end bending moment at A is equal to wrL times the 
area under the curve in Fig. 5 marked 9 = 90, between values of §,/9 
equal to zero and 0.667. The fixed-end bending moment at B is equal to 
wrL times the area between the values of 9,/9 equal to 0.333 and 1.000. . 
Fixed-end twisting moments can then be found from similar curves, and 
those values are combined with the fixed-end bending moments, as 
shown in Eqs. 1, to obtain fixed-end moments about x- and z- axes. 

The writer believes that the use of continuous girders curved in plan 
and the use of continuous frames with curved girders will increase. 
This will be true for bridges and elevated expressways, as well as for 
buildings and foundation frameworks, as witness the new Highbridge 
Interchange at the George Washington Bridge in New York*. The pro- 


‘cedure presented in this paper readily yields a solution for the most 


complicated cases, solution of which is virtually impossible by other 
available methods. 


Correction for Transactions.—In Fig. 5, the horizontal scale should be 
labeled “Values of 


4. “Box Girder on Stilts”, by Howard E. Reid and Phillips H. Lovering, 
Engineering News-Record, Feb. 7, 1952, p. 28. 
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DISCUSSION OF MECHANICS OF MANIFOLD FLOW 
PROCEEDINGS—SEPARATE 258 


JOHN S, McNOWN', M. ASCE.— Most significant among the numer- 
ous contributions of the discussers is the convincing evidence of the 
broad interest and wide range of applicability for studies of manifold 
flow. Applications specifically illustrated were the rapid sand filter 
mentioned by Mr. Barton, the sprinkler irrigation system described in 
useful detail by Messrs. Carstens and Roberson, and the lock manifold 
systems discussed by Mr. Hartigan. Interest is also demonstrated by 
the number of experimenters who have studied the problem. To be 
added to the list referred to in the paper are the Japanese experiments 
mentioned by Mr. Powell, and the Stanford work described by Messrs. 
Vennard and Dentoni. 

One of the most troublesome problems in any experiment is the 
assessment of the probable accuracy of the data. Mr. Barton’s fuller 
explanation of the procedures followed and the nature of the results is 
therefore welcome. In particular, he clearly expressed the importance 
of determining independently for the various pipe sections the piezo- 
metric gradient as a function of the discharge. This measure undoubt- 
edly reduced both the scatter of the data and the absolute error of the 
trends established. 

Additional results and alternative methods of presenting the data 
were contributed by Messrs. Escobar and Niaz, who also had partici- 
pated in the experiments on which the paper was based. Mr. Escobar 
demonstrated convincingly the way in which the transverse velocity dis- 
tribution causes the apparently negative head losses of Fig. 4, and ob- 
served the profile of the contracted jet in the lateral. In Fig. 4 of his 
discussion he illustrated the good correspondence between experiment, 
empiricism, and theory. Mr. Niaz presented in his discussion alterna- 
tive plots, Figs. A and B, for the combining flow. These show (1) the 
degree of similarity of the loss to that in an abrupt contraction, and (2) 
the efficacy of plotting the head loss for the lateral flow against the 
ratio. of the final to the initial velocity for this flow. Although the lat- 
ter plot introduces a degree of system to the results, the data in Fig. B 
appear on close inspection to be three curves, each involving a transi- 
tion from the lower envelope curve to the upper one. The results for 
this type of flow are evidently difficult to depict in any completely satis- 
factory manner. 

The question of the minimum length of conduit or lateral to establish 
the intended uniform flow was raised in two of the discussions. The 
velocity distribution approaches its ultimate form asymptotically with 


1. Prof., Univ. of Michigan, Dept. of Eng. Mechanics, Ann Arbor, Mich.; 
formerly Associate Director, Iowa Inst. of Hydr. Research, State Univ. 
of Iowa, Iowa City, Iowa. 
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distance downstream from the point of separation or major disturbance. 
Selection of the distance required for the branch point no longer to have 
an effect depends very much upon the purpose for which the determina- 
tion is to be made. 

The results shown by Mr. Hartigan are obviously, as he states, sub- 
ject to action related as much to the jet velocity as to the mean velocity 
for the section, the lateral lengths involved being 2 or 3 lateral diame- 
ters. In a distance of 5 to 10 diameters the primary characteristics of 
the separation and high-velocity jet disappear. The relevant questions 
asked by Messrs. Carstens and Roberson can be given only a provi- 
sional answer. For a sprinkler unit to operate in a normal manner, 
again a distance downstream of only a few diameters would be necessary. 
Aside from the indication of the major flow conditions, the proposal to 
use the pressure taps as an indication of the degree to which the flow 
has been established is not reliable. The performance of a subsequent 
lateral, which was tested thoroughly by Yanes (4), was surprisingly de- 
pendent upon small changes in the velocity distribution. These changes 
cannot be detected accurately without a large number of carefully made 
piezometers, and, even then, the significance of a small change in pres- 
sure would probably have been underestimated without the direct meas- 
urements of lateral performance. 

Detection by Messrs. Vennard and Dentoni of the error both in label- 
ling and in discussing the comparison of the Munich and Iowa experi- 
ments in Fig. 4 is greatly appreciated. Although most of the results 
correspond quite well, those for this case certainly do not. Tre differ- 
ences are not large, but the trends are opposite. The comparable data 
from Stanford, generally confirming the Iowa test results, are thus par- 
ticularly valuable; the summary plots presented by Messrs. Vennard 
and Dentoni should also be useful. There are several minor questions 
about these, particularly Fig. 4. The curves should logically approach 
the limiting value of unity as Q,/Q approaches zero. The curve for 
D,/ D = 1.0 is too high at the other end of the scale (Q,/Q = 1.0). The 
suggested value of 1.5 is considerably greater than either the Munich 
value 1.3 or the Iowa value of 1.15. No data were given for D,/ D=1in 
the Stanford tests. If this curve were lowered as it apparently should 
be, the indicated trend of the curves would be improved. Omission of 
numerical values at either end of the scales for the various figures is 
needlessly confusing to the reader, particularly for the logarithmic 
scales. 

Also welcome is the corroboration provided by the additional data of 
the Japanese experimenters for combining flow, referred to in Mr. 
Powell’s discussion. As Mr. Powell pointed out, the correspondence 
among the several sets of data is quite good. For D,/D = 0.5, the Iowa 
and Munich data differ somewhat (Figs. 7 and 8), and in these instances 
the data from Japan are very close to those from Iowa. 

No matter how much care is given to the averaging of the various 
sets of data, the result will not be a final answer to all problems of 
manifold flow. In the paper it was emphasized that the data presented 
are for a single geometrical form of branch point or junction, and that 
the results obtained are representative rather than general. The variety 
of results presented by Mr. Hartigan for various lock ports well 
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lllustrates the variability with port geometry of one characteristic of 
the flow. For this reason, extreme refinement either in the processing 
of the data or in its application to any such design problem as the 
sprinkler system is unjustified. 

Various writers attach importance to different phases of symbology. 
Although the word friction is a misnomer for describing any kind of 
head loss in fluid flow, the symbol hy for head loss is well known, and a 
lower-case subscript is more suitable to the typewriter than is a capi- 
tal one. The terminology involving multiple subscripts, which Messrs. 
Vennard and Dentoni propose, seems needlessly complex. The system 
employed in the paper has the advantages of simplicity and of a maxi- 
mum amount of carry-over between the converging and diverging types 
of flow. 

By way of conclusion, the writer is gratified by the interest shown 
in discussions equivalent in length to nearly twice that of the paper it- 
self, and by the fact that three of those who participated in the Iowa 
tests contributed additional findings from their thesis investigations. 
The goals of the Society in fostering thorough discussion were met ad- 
mirably by the discussers in that varied points of view were expressed, 
additional results were presented, and various aspects of the findings 
were given special emphasis. The representative results obtained, re- 
inforced in the discussions, should find application to a wide variety of 
practical problems. 
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DISCUSSION OF INELASTIC BEHAVIOR OF REINFORCED 
CONCRETE MEMBERS SUBJECTED TO 
SHORT TIME STATIC LOADS 
PROCEEDINGS—SE PARATE 286 


L. H. N. LEE’, A. M. ASCE.—The author wishes to thank the dis- 
cussers for their interests in the paper. Mr. Bresler’s view that the 
precise shape of the flexural stress-strain curve of concrete is pri- 
marily of academic interest is correct as the variations of properties 
of concrete are large. However, there is no doubt that a certain stress- 
strain relationship, which can be mathematically expressed and manip- 
ulated, representing the average properties of concrete is needed. For 
studying the behavior of statically indeterminate structures, buckling of 
columns and other problems, the strength-deformation relationships are 
required; hence, in turn, the stress-strain relationship should be es- 
tablished. 

It is confirmed that the stress-strain curve for concrete tested in 
standard compression cylinder test is a parabolic curve. It is known 
that the state of stress at the compression side of a beam is identical 
with the state of stress at a point inside a compression cylinder. The 
stress-strain relationships of two elements made of identical material 
and subjected to the same state of stress for the same duration obvious- 
ly can not be different. It has also beun early observed by A. N. Talbot 
as he stated in 1906:(16) «, . . . Various curves have been proposed to 
represent the stress-deformation relation for concrete, but the parabola 
is the most satisfactory general representation. Frequently, the parab- 
ola expresses the relation almost exactly and in nearly every case the 
parabolic relation will fit the stress-deformation diagrams very closely 
through the part which is ordinarily developed in beams... .” 

The author has compared the deflections of eccentrically loaded col- 
umns tested by Hognestad, (21) with the theoretical deflections computed 
by the Equations (15) to (29). A total of 95 square and cylindrical col- 
umns were used in the comparison. The theoretical and actual deflec- 
tions agree reasonably well as shown in Fig. 16. 

The reinforcing bars used in the author’s experiments were made of 
mild steel having a sharp yield point. This type of steel is commonly 
used in construction. The stress-strain curve of the steel is almost 
exactly represented by the idealized stress-strain curve shown in Fig. 
4. In this case, the proportional limit and the yield point are very close 
to each other. 

Mr. Bresler’s example of using Equations (3) and (4) is not correct. 
Equation (3) is for the range of 0<€.< and Equation (4) is for the 
range of €y< €s<€yp. The point€, =€y is a point of discontinuity in 


I. Asst. Prof., Dept. of Eng. Mechanics, Univ. of Notre Dame, Notre 
Dame, Ind. 


16. Refer to footnote on page 286-10. 
21. Refer to footnote on page 286-13. 
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first derivative. The values of de,/d€, for the left side and for the 
right side of the yield point are not the same. 

Equation (3') suggested by Mr. Bresler is useful for the range when 
steel stresses are between proportional limit and yield point. 

It is correct that calculation of flexural stress-strain relation in- 
volves numerical differentiation which may be very sensitive to the ac- 
curacy of experimental measurements as pointed out by Mr. Bresler. 

It has been experienced by the author that tests should be carried out 
with uniform increment of load in uniform time interval. Total time of 
testing should be short and special attention is needed at the ultimate 
load stage. It has also been the author’s experience that use of a simply- 
supported beam will not give a stress-strain relation over a large range 
of strain, because the rapid collapse of the beam after it has absorbed a 
maximum of strain energy prevents the taking of an adequate number of 
accurate readings. 

It should be pointed out that the test results with the strains less than 
6 x 10™* shown in Fig. 15 as presented by Mr. Bresler are not valid for 
the assumption of no tensile stress in concrete is not valid in that range. 

The method of obtaining flexural stress-strain relationship suggested 
by Mr. Ashdown is questionable; because under the suggested condition, 
there is no way to preserve the condition of plane remains plane which 
is a basic assumption in forming the stress-strain relationship. 

The readings shown in Table I stopped at 98% of the yielding moment 
of the.section and 844% of the collapse load of the restrained beam as 
clearly shown in Fig. 13. The strains beyond that load were not ob- 
tained, because the strains changed swiftly at that stage. The distribu- 
tion of stresses as estimated by Mr. Ashdown is actually shown in 
Fig. 6. 


The author does not agree with Mr. Ashdown that the ordinary engi- 
neer will be frightened by the complicated formulae given in the paper. 
However, if Mr. Ashdown’s estimate of the “Ordinary Engineer” is val- 
id, there remains the possibility of putting the results in tabular and 
graphical form, convenient for computation. 
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DISCUSSION OF IMPULSIVE MOTION 
OF ELASTO-PLASTIC BEAMS 


H. H. BLEICH’, M. ASCE and M. G. SALVADORI’, M. ASCE.—While 
the discusser concedes in his opening paragraph that the writers “have 
made a valuable contribution” in their paper, he objects later on to the 
results contained in it, which he calls variously “certainly wrong,” “er- 
roneous” and “untenable”. The reasoning on which these objections are 
based is not intended, apparently, to attack the theoretical part of the 
paper, but rather the manner of its application to the particular numeri- 
cal examples presented. The objections are based on differences be- 
tween the writers’ results and those of a so-called “correct rigid- 
plastic” analysis. The writers will show in this closure that this latter 
analysis is not applicable to their examples and that, therefore, the dis- 
cusser has no case. 

It is the writers’ contention that, for reasons set out below, the man- 
ner of handling the examples is justified by the physical facts, and it will 
be shown here that the results of the so-called “correct rigid-plastic” 
analysis do not approximate the behavior of the actual elasto-plastic 
beam, as claimed by the discusser. This claim will be proved unjusti- 
fied because the energy criterion used by the discusser is incomplete 
and a second essential criterion must be taken into consideration to de- 
cide whether a rigid-plastic solution of an initial velocity problem is a 
realistic approximation to an elasto-plastic case. 

In order to disprove the discusser’s claim that the so-called “correct 
rigid-plastic” solution is a realistic approximation for an actual elasto- 
plastic beam, the nature of the response of the beam for small values of 
the time t in both cases will be considered. As Miss Cotter’s paper (to 
which the discusser refers) is unpublished and unavailable at the time 
of writing and as the discusser does not state otherwise, it is assumed 
that Miss Cotter’s approach is essentially the one used by the discusser 
in various papers already published, but adapted to the case of an initial 
velocity distribution without external forces. This means that a bar sub- 
jected to an initial velocity distribution Yo(x) will yield at t = 0 in all 


locations x where the curvature of the bar ae) =y'(x) changes 


with time, i.e. where 7" (x) 4 0; and that in all locations where yield oc- 
curs the bending moment is equal to the capacity moment M,. However, 
the sign of the moment, and the sign of the permanent curvature develop- 
ing in any location depend on the sign of y (x). After the assumed incep- 
tion of yielding at t = 0, but for very small values of t the deformation 
y(x, t) is proportional to y,,(x), or 


y(x,t) =ty,(x) [for very small t) (1) 


1. Prof. of Civ. Eng., Columbia Univ., New York, N. Y. 
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and the time rate of change of the curvature, the sign of which controls 
the direction of the permanent curvature, becomes 


¥"(x) = ¥o"(x) ¢ [for very small t] (2) 


The sign of y"(x) is simply the sign of the second derivative of the ini- 
tial velocity distribution. 

Applying this result to the velocity distribution of the example, Fig. 3 
of the paper, 


Yo 7x 
Vo (x) = (1 + cos 


it is seen that y"(x) is positive for the outer quarters of the beam length, 
but negative in the center half. The so-called “correct rigid-plastic” 
analysis indicates therefore that, for early times, angular changes in 
opposite directions will occur in the various parts of the beam. Ast 
increases and energy is absorbed this behavior will of course be modi- 
fied in a manner which is the subject of Cotter’s paper referred to ear- 
lier. For sufficiently large t the beam will be permanently deformed to 
a shape shown qualitatively in Fig. 1 and will move with a certain final 
velocity y(x) = v,,. The center portion will be bent through a total angle 
6,, while each end will be bent in the opposite direction by angles @,. 
(The angle 9 * used by the discusser is not clearly defined; it may be 
either 6*=6,, or 0*= 6, - 26,. The difference is immaterial in what 
follows.) 

The permanent deformation obtained by the writers is a kink at the 
center, at the assumed yield hinge. This assumption has been made be- 
cause the writers expect the permanent deformation of the beams in the 
examples not to have the character of Fig. 1, but rather of Fig. 2, where 
all permanent curvatures have the same sign, and occur essentially in 
the center portion of the beam. If the deformation is of the character of 
Fig. 2, an assumed yield hinge at the center will clearly give physically 
meaningful values of the total angle @ for the elasto-plastic case. The 
value A. as defined by the writers will be an upper bound for the 
elasto-plastic value; the fact that a rigid-plastic approach was used to 
obtain 6,, does not mean that the result, in general, must be equal to the 
“rigid-plastic” solution of Miss Cotter. For the velocity distribution 
of Eq. (3) the results can not possibly agree because considerations of 
energy require that 6,, = 0, + 26, and therefore that 0 4 6, since 0, 40. 

The crux of the argument is the writers’ contention that the counter- 
flexures obtained in the rigid-plastic approach do not occur in the elasto- 
plastic beams in either of the examples in the paper. Consider the re- 
sponse of the elastic beam (prior to the plastic stage) and the time his- 
tory of the bending moment for small values of t. For sufficiently small 
values of t, the displacements are given again by Eq. (1), and the bend- 
ing moments become 


M(x, t) = - Ely," (x)t [for very small t] (4) 


For the velocity distribution under consideration, Eq. (4) indicates that 

the largest moments, of opposite signs, occur at x = 0 andx =+i. 

However, the points x = + are at the ends of the beam where the mo- 

ments must necessarily vanish; Eq. (4) therefore does not apply, even 

for small values of t if { - Ix) is small. One can, however, visualize the 
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actual bending moments to be the sum of Eq. (4) and of correcting stress 
waves originating at x = + { such that for the sum the condition M(+X, t) 
= 0 remains satisfied.? This suggests that Eq. (4) will apply at points 
near the ends x = + { only if t is smaller than the time required for a 
stress wave to travel the distance { - |x|. In the present case this 
stress wave will decrease all moments near the ends, but increase 
those at the center of the beam. When the time t is larger than i /¢, 

the time required by the waves to travel through the entire beam, Eq. (4) 
becomes inapplicable for all x, and the bending moments will be redis- 
tributed in a manner which need have no similarity to the original curve 
Yo (x). If first yield takes place at a time ty? c, the sign of the bend- 
ing moments at yield and the direction of the permanent curvature will 
no longer be governed by the sign of ¥o" (x). It follows that in the ex- 
amples considered the type of permanent deformation shown in Fig. 1 
can only be obtained if the time to yield is very small compared to Vc. 
If this is not the case the counterflexing moments will have disappeared 
by the time yielding begins. 

The reasoning in the previous paragraph is physically incomplete be- 
cause shear deformations have been neglected; for sufficiently small 
values of the time t this is not permissible. It is not difficult, however 
to go through the same reasoning for the case of the Timoshenko beam 
allowing for shear deformations. Excluding cases where yield in shear 
occurs (which are not contemplated here), one finds again that the yield 
pattern is governed by the sign of y,"(x) only if the time at yield is small 
versus {/c, confirming the conclusion found without shear effects. 

The so-called “correct rigid-plastic” analysis obtains the location 
and direction of yielding from the values of Yo" (x) and one must con- 
clude that this analysis can not be relied upon when the time to yield in 
the elasto-plastic structure is of the order of { /c or larger, because in 
such a case the actual yield pattern is quite different from that obtained 
by the so-called “correct rigid-plastic” analysis. This difference is im- 
portant for the examples used in the paper. (It may not be important for 
other cases, but the question of other cases does not rightly belong here.) 

The examples in the paper may be interpreted as applying to steel 
beams by selecting suitable units, e.g. ft., sec., and kips. The times to 
first yield, 0.2 sec. and 0.016 sec. respectively are not small versus 
\/e ~ 0.01 sec; hence, the so-called “correct rigid-plastic” solution is 
not applicable in either case. It should be clear that the time to yield 
would have to be considerably smaller than } /e to make the rigid- 
plastic solution applicable. 

The reviewer asks also why the writers have not discussed “methods 
for examining all cross-sections of the beam.” No special “methods” 
are required. Once the coordinates q; or Pp; are obtained the moment 
history can be plotted by routine computation, if one desires to do so. 
Some points of such curves were computed, but not included in the paper 
as they did not reveal anything of particular interest. The raising of the 
above question, and the discusser’s emphasis on checking the capacity 
moments of adjoining sections shows a lack of realization of the type of 


2. To be consistent, the beam should be deemed to have rotatory inertia 
and to be infinitely rigid in shear. Eq. (4) does not visualize shear de- 
formations. 
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problem which should and can be treated by normal mode methods in 
general, with or without plasticity. If one desires to treat a problem in 
which information at many points is required with good accuracy, nor- 
mal mode methods are impractical even in the elastic case. On the 
other hand, if one is looking for approximate, yet physically reasonable 
results, mode methods come into their own, because one can restrict the 
computations to a small number of modes. Such a restriction results in 
a loss of detail and in problems of elasto-plastic beams it would be ab- 
surd to pay attention to minor motion of hinges or areas becoming plas- 
tic; any conclusions as to such matters on the basis of a few terms of 
an expression would be quite unrealistic. It is only consistent with the 
limited accuracy obtainable to place a yield hinge in a suitably selected 
strategic position, or in several positions if necessary; one will still 
get useful and meaningful results. This is the reason for the selection 
of a yield hinge at the center only in the examples treated; it represents 
the essential physical effects. If the initial velocity in the examples j 
were increased considerably it might of course be necessary to consider 7 
two additional hinges in the outer quarters of the beam to allow for per- 
manent deformations according to Fig. 1. Cotter’s rigid-plastic solution 
is an excellent basis for the selection of the best location of these hing- 
es. For sufficiently large v, where the rigid-plastic solution is appli- 
cable there is of course no point in even trying to use a raode solution. 

To sum up, it is the writers’ position that both elasto-plastic and 
rigid-plastic methods have their place. The elastc-plastic approach is 
necessary in the very important practical cases where the loads exceed 
the elastic capacity by a modest amount only, since a rigid-plastic anal- 
ysis (which is necessarily simpler) is inapplicable to such cases. j 
Moreover, the writers point out that in cases of initial velocities an ad- . 
ditional criterion may be required for the rigid-plastic solution to be : 
valid. (In the original publication by Lee and Symonds the energy cri- 
terion was not really proposed for cases of initial velocity, since these j 
cases were not treated there at all.) eg 

In the very last paragraph of the discussion the discusser raises the i 
question of simple beams on which subject he is preparing a paper. It 
seems best to postpone the consideration of this question until this paper 
is available. 


Postscript 


After submission of the preceding closure, an advance copy of Miss 
Cotter’s paper to which the discusser refers was made available to the 
writers. One of the results of this paper, not mentioned by the discus- 
ser, is of considerable interest to this closure and is therefore added 
here. 

Miss Cotter found that the deformation of the bar is of the type shown 
in Fig. 1 but that the total permanent deformation 6, near the center of 
the bar consists of two parts: one part is a kink at the center, the other 
a bend distributed over a length of about //6, 

In comparing Miss Cotter’s results and those of the writers, the dis- 
cusser has quoted only the value of the kink-angle at the center; this 
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angle differs from the kink-angle found by the writers on the basis of 
the single hinge analysis by the stated amounts. Miss Cotter’s paper 
also contains the values for the distributed bending in the vicinity of the 
center, and the total angles found by Miss Cotter and by the writers 
differ only very slightly: 


Initial velocity v, = 11.1 67.4 
Miss Cotter: 2.64 Vo 10-5 1.60 Yo 10-2 
Bleich and Salvadori: 2.67 Vo 10-5 1.61 Yo 10-? 


The “so-called correct rigid-plastic” analysis differs therefore 
from the writers’ upper limit solution essentially only by the fact that 
some of the deformation is distributed over a short distance rather than 
concentrated as a kink. It is worth recalling that the occurrence of any 
kink is only due to the simplification in the moment-curvature relation- 
ship, and that for purposes of application to any actual structure only 
the total angle is physically meaningful. This indicates that the differ- 
ences between the two rigid-plastic analyses are pretty close to trivial. 

The above comparison is not intended to withdraw the writers’ argu- 
ment in the main body of this closure, stressing that the “so-called 
correct rigid-plastic” analysis is not applicable ‘in either of the ex- 
amples because the additional criterion stated is not satisified. 

In conclusion, the writers believe to have shown that the discusser’s 
criticism of the results of their paper is not justified. 


CORRECTIONS for Transactions.—On the front cover page, the initials 
of M. G. Salvadori have been printed incorrectly. In Eq. 13, the double- 
prime superscript should be outside the closing parenthesis. In Fig. 8, 


the expression at the top of the ordinate scale should be _ 10°60 
Yo 


=a 
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